Dynamics of Symmetry Breaking and Tachyonic Preheating 
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We reconsider the old problem of the dynamics of spontaneous symmetry breaking using 3d lat- 
tice simulations, and develop a theory of tachyonic preheating, which occurs due to the spinodal 
instability of the scalar field. Tachyonic preheating is so efficient that symmetry breaking typically 
completes within a single oscillation of the field distribution as it rolls towards the minimum of 
its effective potential. As an application of this theory we consider preheating in the hybrid infla- 
tion scenario, including SUSY-motivated F-term and D-term inflationary models. We show that 
preheating in hybrid inflation is typically tachyonic and the stage of oscillations of a homogeneous 
component of the scalar fields driving inflation ends after a single oscillation. Our results may also 
be relevant for the theory of the formation of disoriented chiral condensates in heavy ion collisions. 
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I. INTRODUCTION 

Spontaneous symmetry breaking is a basic feature of 
all realistic theories of elementary particles. In the sim- 
plest models, this instability appears because of the pres- 
ence of tachyonic mass terms such as —m 2 (f> 2 /2 in the 
effective potential. As a result, long wavelength quan- 
tum fluctuations <f>k of the field (f> with momenta k < m 
grow exponentially, tfik ~ exp(iV m 2 — k 2 ), which leads 
to spontaneous symmetry breaking. 

This process may occur gradually, as in the theory 
of second order phase transitions, when the parameter 
m 2 slowly changes from positive to negative and the de- 
ree of symmetry breaking gradually increases in time 
Sometimes the symmetry breaking occurs discontin- 
uously, due to a first order phase transition Q . But there 
is also another possibility, which we will study in this pa- 
per: The tachyonic mass term may appear suddenly, on 
a time scale that is much shorter than the time required 
for symmetry breaking to occur. This may happen, for 
example, when the hot plasma created by heavy ion col- 
lisions in the 'little Big Bang' suddenly cools down ||. 
A more important application from the point of view of 
cosmology is the process of preheating in the hybrid in- 
flation scenario , where inflation ends in a 'waterfall' 
regime triggered by tachyonic instability. 

During the last few years we have learned that particle 
production by an oscillating scalar field may occur within 
a dozen oscillations due to the nonperturbative process 
of particle production called preheating || . Usually pre- 
heating is associated with broad parametric resonance in 
the presence of a coherently oscillating inflaton field ||, 
but other mechanisms are also possible, see e.g. 

In this paper we will concentrate on what we call tachy- 



onic preheating, which occurs due to tachyonic (spin- 
odal) instabilities in the field responsible for the sym- 
metry breaking. The process of symmetry breaking has 
been studied before by advanced methods of perturba- 
tion theory, see e.g. || and references therein. However, 
spontaneous symmetry breaking is a strongly nonlinear 
and nonperturbative effect. It usually leads to the pro- 
duction of particles with large occupation numbers in- 
versely proportional to the coupling constants. As a re- 
sult the perturbative description, including the Hartree 
and 1/N approximations, has limited applicability. It 
does not properly describe rescattering of created parti- 
cles and other important features such as production of 
topological defects. 

For further theoretical understanding of the issue one 
should go beyond perturbation theory. Fortunately, dur- 
ing the last few years new methods of lattice simulations 
have been developed. They are based on the observation 
that quantum states of bose fields with large occupa- 
tion numbers can be interpreted as classical waves and 
their dynamics can be fully analyzed by solving relativis- 
tic wave equations on a lattice p0| , pl| . In our paper we 
will further develop these methods, including effects of 
renormalization, and apply them to the investigation of 
spontaneous symmetry breaking and tachyonic preheat- 
ing. A significant advantage of these methods as com- 
pared to other lattice simulations of quantum processes 
is that the semi-classical nature of the effects under inves- 
tigation allows us to have a clear visual picture of all the 
processes involved. That is why this paper is accompa- 
nied by several computer generated movies that show the 
development of symmetry breaking in various models. 

We will show that tachyonic preheating can be ex- 
tremely efficient. In many models it leads to the transfer 
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of the initial potential energy density V(0) into the en- 
ergy of scalar particles within a single oscillation. Con- 
trary to standard expectations, the first stage of preheat- 
ing in hybrid inflation is typically tachyonic, which means 
that the stage of oscillations of a homogeneous compo- 
nent of the scalar fields driving inflation either does not 
exist at all or ends after a single oscillation. A detailed 
description of our results will be given elsewhere . 

II. TACHYONIC INSTABILITY AND 
SYMMETRY BREAKING 

Symmetry breaking occurs due to tachyonic instability 
and may be accompanied by the formation of topological 
defects. Here we will consider two toy models that are 
prototypes for many interesting applications, including 
symmetry breaking in hybrid inflation. 

A. Quadratic potential 

The simplest model of spontaneous symmetry breaking 
is based on the theory with the effective potential 

where A « 1. V(<p) has a minimum at = ±i> and a 
maximum at = with curvature V" = —m 2 . 

The development of tachyonic instability in this model 
depends on the initial conditions. We will assume that 
initially the symmetry is completely restored so that the 
field does not have any homogeneous component, i.e. 
(0) = 0. But then (0) remains zero at all later stages, 
and for the investigation of spontaneous symmetry break- 
ing one needs to find the spatial distribution of the field 
0(x, t). To avoid this complication, many authors assume 
that there is a small but finite initial homogeneous back- 
ground field 0(i) , and even smaller quantum fluctuations 
6(j>(x,t) that grow on top of it. This approximation may 
provide some interesting information, but quite often it 
is inadequate. In particular, it does not describe the cre- 
ation of topological defects, which, as we will sec, is not a 
small nonperturbative correction but an important part 
of the problem. 

For definiteness, we suppose that in the symmetric 
phase = there are usual quantum fluctuations of the 
massless field with the mode functions _i e - lfc *+ lfca: anc [ 

then at t = we 'turn on' the term — to 2 2 /2 correspond- 
ing to the negative mass squared — to 2 . The modes with 
k = \k\ < to grow exponentially so the dispersion of these 
fluctuations can be estimated as 

m 

l 6( f?) = JL / dkke 2t ^^ . (2) 
\n l J 
o 



To get a qualitative understanding of the process of 
spontaneous symmetry breaking, instead of many grow- 
ing waves with momenta k < to in (||) let us consider first 
a single sinusoidal wave 50 = A(t) cos kx with k ~ m and 
with initial amplitude ~ S in one-dimensional space. 
The amplitude of this wave grows exponentially until its 
becomes 0{v) ~ m/VX- This leads to the division of the 
universe into domains of size C(m _1 ) in which the field 
changes from 0(v) to O(-v). The gradient energy den- 
sity of domain walls separating areas with positive and 
negative will be ~ fc 2 (50 2 = 0(m A /X). This energy is 
of the same order as the total initial potential energy of 
the field V(0) = m 4 /4:X. This is one of the reasons why 
any approximation based on perturbation theory and ig- 
noring topological defect production cannot give a cor- 
rect description of the process of spontaneous symmetry 
breaking. 

Thus a substantial part of the false vacuum energy 
V(0) is transferred to the gradient energy of the field <p 
when it rolls down to the minimum of V (0). Because the 
initial state contains many quantum fluctuations with 
different phases growing at a different rate, the resulting 
field distribution is very complicated, so it cannot give 
all of its gradient energy back and return to its initial 
state = 0. This is one of the reasons why spontaneous 
symmetry breaking and the main stage of preheating in 
this model may occur within a single oscillation of the 
field 0. 

Meanwhile if one were to make the usual assumption 
that initially there exists a small homogeneous back- 
ground field <C v with an amplitude greater than the 
amplitude of the growing quantum fluctuations (50, so 
that m/27r <C < m/VX, one would find out that when 
falls to the minimum of the effective potential the gra- 
dient energy of the fluctuations remains relatively small. 
One would thus come to the standard conclusion that 
the field should experience many fluctuations before it 
relaxes near the minimum of V(0). To avoid this er- 
ror, we need to perform a complete study of the growth 
of all tachyonic modes and their subsequent interaction 
without making this simplifying assumption about the 
existence of the homogeneous field 0. 

Consider the tachyonic growth of all fluctuations with 
k < to, i.e. those that contribute to (<50 2 ) in Eq. (|^). 
This growth continues until \J (<50 2 ) reaches the value 
~ v/2, since at ~ the curvature of the effective 

potential vanishes and instead of tachyonic growth one 
has the usual oscillations of all the modes. This happens 
within the time f * ~ In ^- . The exponential growth 
of fluctuations up to that moment can be interpreted as 
the growth of the occupation number of particles with 
k <C to. These occupation numbers at the time i* grow 
up to 

/ 7T 2 \ 7T 2 

7ifc ~ exp(2TOi„) ~ exp I In — I = — 3> 1 ■ (3) 
One can easily verify that t* depends only logarithmi- 
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cally on the choice of the initial distribution of quantum 
fluctuations. For small A the fluctuations with k <C m 
have very large occupation numbers, and therefore they 
can be interpreted as classical waves of the field 0. 

The dominant contribution to (6(f) 2 ) in Eq. (^) at 
the moment i* is given by the modes with wavelength 
h ~ 271-fc- 1 ~ V^Trm-Mn 172 (Cn 2 /X) > m" 1 , where 
C = 0(1). As a result, at the moment when the fluc- 
tuations of the field <f> reach the minimum of the effective 
potential, v/ (0 2 ) ~ v, the field distribution looks rather 
homogeneous on a scale I < /*. On average, one still 
has ((f>) = 0. This implies that the universe becomes 
divided into domains with two different types of spon- 
taneous symmetry breaking, cj> ~ ±v. The typical size 
of each domain is Z»/2 ~ -3= m^ 1 In 1 / 2 ^j- 7 which differs 

only logarithmically from our previous estimate m . At 
later stages the domains grow in size and percolate (eat 
each other up), and spontaneous symmetry breaking be- 
comes established on a macroscopic scale. 

Of course, these are just simple estimates which should 
be followed by a detailed quantitative investigation. 
When the field rolls down to the minimum of its effective 
potential, its fluctuations scatter off each other as clas- 
sical waves. It is difficult to study this process analyti- 
cally, but fortunately one can do it numerically using the 
method of lattice simulations developed in Jl^,|ll[. We 
performed our simulations on lattices with either 128 3 
and 256 3 gridpoints. A detailed description of our calcu- 
lations will be given in |l2j ; here we will only present our 
main results. 

Figure [l] illustrates the dynamics of symmetry breaking 
in the model ([!]). It shows the probability distribution 
P((j),t), which is the fraction of the volume containing 
the field ^ at a time t if at t = one begins with the 
probability distribution concentrated near cj> = 0, with 
the quantum mechanical dispersion (^). 

As we see from this figure, after the first oscillation the 
probability distribution P(<fi, t) becomes narrowly con- 
centrated near the two minima of the effective potential 
corresponding to (f> = ±v. In this sense one can say that 
symmetry breaking completes within one oscillation. To 
demonstrate that this is not a strong coupling effect, we 
show the results for the model (|j) with A = 1CP 4 . Note 
that only when the distribution stabilizes and the do- 
mains become large can one use the standard language 
of perturbation theory describing scalar particles as ex- 
citations on a (locally) homogeneous background. That 
is why the use of the nonperturbative approach based on 
lattice simulations was so important for our investigation. 

The growth of fluctuations in this model is shown 
in Fig. [2| It shows how fluctuations grow in a two- 
dimensional slice of 3D space. Maxima correspond to 
domains with <p > 0, minima correspond to domains with 
6 < 0. 
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FIG. 1. The process of symmetry breaking in the model 
(^) for A = 10~ 4 . In the beginning the distribution is very 
narrow. Then it spreads out and shows two maxima which 
oscillate about <f> — ±v with an amplitude much smaller than 
v. These maxima never come close to the initial point <f> = 0. 
The values of the field are shown in units of v. 




FIG. 2. Growth of quantum fluctuations in the process of 
symmetry breaking in the quadratic model (Q). 

The dynamics of spontaneous symmetry breaking 
in this model is even better illustrated by the com- 
puter generated movie is included in this submis- 
sion as a file l.gif, which can be found also at 



http://physics.stanford.cdu/gfcldcr/hybrid/l.gil. It con- 
sists of an animated sequence of images similar to the one 
shown in Fig. 0. These images show the whole process 



3 



of spontaneous symmetry breaking from the growth of 
small gaussian fluctuations of the field <f> to the creation 
of domains with <j> = ±v. 

Similar results are valid for the theory of a multi- 
component scalar field (f>i with the potential ([!]). For 
example, the behavior of the probability distribution 
P(4>i,(j)2,t) in the theory of a complex scalar field 
cf> = (4>i + i02)/\/2 is shown in Fig. [|. As we 
see, after a single oscillation this probability distribu- 
tion has stabilized at \(f>\ ~ v. A computer gener- 
ated movie illustrating this process is included in this 
submission as a file 2.gif, which can be found also at 
http://physics.stanford.cdu/gfelder/hybrid/2.gii. 
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FIG. 3. The process of symmetry breaking in the model 
(0) for a complex field (f>. The field distribution falls down to 
the minimum of the effective potential at \cj>\ = v and expe- 
riences only small oscillations with rapidly decreasing ampli- 
tude \A(j>\ < v. 



B. Cubic potential 

Another important example of tachyonic preheating is 
provided by the theory 
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(4) 



This potential is a prototype of the potential that appears 
in descriptions of symmetry breaking in F-term hybrid 
inflation 

The first question to address concerns the initial am- 
plitude of the tachyonic modes in this model. This is 
nontrivial because to 2 (</>) = —2\v<p + 3A</> 2 vanishes at 
4> = 0. However, eq. (||) implies that scalar field fluc- 
tuations with momentum ~ k have initial amplitude 
A - - They enter a self-sustained tachyonic 



_k 
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The average initial amplitude of the growing tachyonic 
fluctuations with momenta smaller than is 

These fluctuations grow until the amplitude of 5<p be- 
comes comparable to 2v/3, and the effective tachyonic 
mass vanishes. At that moment the field can be repre- 
sented as a collection of waves with dispersion \J {8(f) 2 ) ~ 
v, corresponding to coherent states of scalar particles 

with occupation numbers nu ~ ("T - ) 

Because of the nonlinear dependence of the tachy- 
onic mass on a detailed description of this process 
is more involved than in the theory ([!]). Indeed, even 
though the typical amplitude of the growing fluctuations 
is given by (||), the speed of the growth of the fluctu- 
ations increases considerably if the initial amplitude is 
somewhat bigger than (||). Thus even though the fluc- 
tuations with amplitude a few times greater than (|^) 
are exponentially suppressed, they grow faster and may 
therefore have greater impact on the process than the 
fluctuations with amplitude (^|). Low probability fluc- 
tuations with 5<j) 3> <5</>rms correspond to peaks of the 
initial Gaussian distribution of the fluctuations of the 
field <p. Such peaks tend to be spherically symmetric 
|]l5| . As a result, the whole process looks not like a uni- 
form growth of all modes, but more like bubble produc- 
tion (even though there are no instantons in this model). 
The results of our lattice simulations for this model are 
shown in Fig. |[ The bubbles (high peaks of the field 
distribution) grow, change shape, and interact with each 
other, rapidly dissipating the vacuum energy V(0). A 
computer generated movie illustrating this process is in- 
cluded in this submission as a file 3.gif. It can be found 
also at http://physics.stanford.cdu/gfelder/hybrid/3.gii. 
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regime if k 2 < \m 2 cS \ = 2Xvy/{S^) - i.e. if k < |f . 



FIG. 4. Fast growth of the peaks of the distribution of the 
field d> in the cubic model (Q). It should be compared with 
Fig. Q for the quadratic model (m) . 
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Fig. @ shows the probability distribution P{tf>, t) in the 
model (Q). As we see, in this model the field also relaxes 
near the minimum of the effective potential after a single 
oscillation. 
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FIG. 5. Histograms describing the process of symmetry 
breaking in the model (^) for A = 1CP 2 . After a single os- 
cillation the distribution acquires the form shown in the last 
frame and after that it practically does not oscillate. 

One should note that numerical investigation of this 
model involved specific complications due to the neces- 
sity of performing renormalization. Lattice simulations 
involve the study of modes with large momenta that are 
limited by the inverse lattice spacing. These modes give 
an additional contribution to the effective parameters of 
the model. In the simple model (|l|) these corrections were 
relatively small, but in the cubic model they induce an 
additional linear term Av0(0 2 ). This term should be sub- 
tracted by the proper renormalization procedure, which 
brings the effective potential back to its form (Q). 

For completeness we would like to mention also that in 
the theory with the quartic potential V = V(0) — \\(j> A 
the decay of the symmetric phase occurs via tunneling 
and the formation of bubbles, even though there is no 
barrier between = and 0=^0 |1(J. In this case 
quantum tunneling can be heuristically interpreted as a 
building up of stochastic fluctuations 5cf> In this 

respect the character of tachyonic instability for the cu- 
bic potential is intermediate between the quadratic and 
quartic potentials. 



III. TACHYONIC PREHEATING IN HYBRID 
INFLATION 

The results obtained in the previous section have im- 
portant implications for the theory of reheating in the 
hybrid inflation scenario. The basic form of the effective 
potential in this scenario is B 



V(cf>,a) = ~(a 2 



b 2 a 2 + ^m 2 (, 



(6) 



The point where cj> — <j) c — M/g and a — is a bi- 
furcation point. Here M = \f\v. The global mini- 
mum is located at = and |er| = v. However, for 
(f) > <fi c the squares of the effective masses of both fields 
m 2 = g 2 (f> 2 — \v 2 + 3\<7 2 and mi = m 2 +g 2 <j 2 are positive 
and the potential has a valley at a = 0. Inflation in this 
model occurs while the 4> Reld rolls slowly in this valley 
towards the bifurcation point. When reaches C , infla- 
tion ends and the fields rapidly roll towards the global 
minimum at = 0, |er| = v. If a is a real one-component 
scalar, this may lead to the formation of domain walls. 
To avoid this problem, we assume that a is a complex 
field. In this case symmetry breaking after inflation pro- 
duces cosmic strings instead of domain walls Q . 

In realistic versions of this model the mass m is ex- 
tremely small, as well as its initial velocity 0. The fields 
fall down along a certain trajectory 0(t),cr(i) in such a 
way that initially this trajectory is absolutely flat, then 
it rapidly falls down, and then it becomes flat again near 
the minimum of V(0, cr). This implies that the curvature 
of the effective potential along this curve is initially neg- 
ative. Therefore the fields should experience tachyonic 
instability along the way. 

The decay of the homogeneous inflaton field and pre- 
heating in hybrid inflation were considered in two papers: 
for the simplest non-supersymmetric scenario with a va- 
riety of parameters j^] and for a SUSY F-term model . 
Both papers were focused on the possibility of parametric 
resonance. However, in [^) it was also pointed out that 
for g 2 ^> A the field a falls down only when the field 
reaches some point ^ < As a result, the motion of 
the field a occurs just like the motion of the field in the 
theory ([!]). In this case one has a tachyonic instability 
and the fields relax near the minimum of V((f>, a) within 
a single oscillation For all other relations between 
g 2 and A the fields follow more complicated trajectories. 
One might expect that the fields would in general expe- 
rience many oscillations, which might or might not lead 
to parametric resonance HQ- 

We performed an investigation of preheating in hybrid 
inflation in the model (^|) with two scalar fields (one real 
and one complex) and in SUSY-motivated F-term and 
D-term inflation models with three complex fields. We 
used methods similar to those that we applied in the pre- 
vious section to the investigation of spontaneous symme- 
try breaking, including 3D lattice simulations. We have 
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found that efficient tachyonic preheating is a generic fea- 
ture of the hybrid inflation scenario, which means that 
the stage of oscillations of the quasi-homogeneous com- 
ponents of the scalar fields driving inflation is typically 
terminated by the backreaction of fluctuations. Here we 
report the qualitative results of our findings, leaving the 
technical details for the longer paper |12| . 

In particular, Fig. |^ shows the process of spontaneous 
symmetry breaking in the theory (Q) for g 2 = I0~ 4 , 
A = 1(T 2 , M = 10 15 GeV. The probability distribu- 



tion oscillates along the ellipse g 
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As before, it relaxes near the minimum of the effec- 
tive potential within a single oscillation. A computer 
generated movie illustrating this process is included in 
this submission as a file 4.gif. It can be found also at 



http://physics. Stanford. edu/gfelder/hybrid/4.gil. 
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FIG. 6. The process of symmetry breaking in the hybrid 
inflation model (H) for g 2 <C A. The field distribution moves 
along the ellipse g (f> 2 + Xa 2 — g 2 (f>c from the bifurcation point 
<b = d>c, a = 0. 



The theory of preheating in D-term inflation for vari- 
ous relations between g 2 and A [[l8) is very similar to the 
theory discussed above. Meanwhile, in the case g 2 = 2A 
the effective potential (^|) has the same features as the ef- 
fective potential of SUSY- inspired F-term inflation [ p~3| . 
In this scenario the fields <f> and ct fall down along a sim- 



ple linear trajectory |I4|, so that instead of following each 
of these fields one may consider a linear combination of 
them and find the effective potential in this direction. 
This effective potential has exactly the same shape as 
our cubic potential (0). Thus all of the results which 
we obtained for tachyonic preheating in the theory (Q) 
should be valid for F-term inflation as well, with minor 
modifications due to the presence of additional degrees of 
freedom that can be excited during preheating. Indeed, 
we were able to confirm these conclusions by lattice sim- 



ulations of the F-term and D-term models. 

Thus we see that tachyonic preheating is a typical fea- 
ture of hybrid inflation. The production of bosons in this 
regime is nonperturbative, very fast, and efficient, but it 
is usually not related to parametric resonance. Instead 
it is related to the production and scattering of classical 
waves of the scalar fields. Of course, one should keep 
in mind that there may exist some particular versions of 
hybrid inflation in which tachyonic preheating is ineffi- 
cient, e.g. because of fast motion of the field <j> near the 
bifurcation point. 

So far we have discussed tachyonic preheating in the 
inflaton sector of hybrid models, which leads to the decay 
of the homogeneous fields and the excitation of their fluc- 
tuations. Preheating in the non- inflaton sector and the 
subsequent development of equilibrium in hybrid models 
were considered in |fl9|| . Light bosonic fields interacting 
with scalars from the inflaton sector are dragged into the 
process of preheating. Excitations of these fields rapidly 
acquire large occupation numbers and further evolve into 
equilibrium together with the scalars from the inflaton 
sector. 

The tachyonic nature of preheating in hybrid inflation 
implies, in particular, that instead of the production of 
gravitinos by a coherently oscillating field |^,^f| , in hy- 
brid inflation one should study gravitino production due 
to the scattering of classical waves of the scalar fields pro- 
duced by tachyonic preheating. Our results may also be 
important for the theory of the generation of the baryon 
asymmetry of the universe at the electroweak scale |22j . 

From a more general point of view, however, the most 
important application of our results is to the general the- 
ory of spontaneous symmetry breaking. This theory con- 
stitutes the basis of all models of weak, strong and elec- 
tromagnetic interactions. The new methods developed 
during the last few years in application to the theory of 
reheating after inflation have been applied in this paper 
to the theory of spontaneous symmetry breaking. These 
methods have for the first time allowed us not only to 
calculate correlation functions and spectra of produced 
particles, but to actually see the process of spontaneous 
symmetry breaking and to reveal some of its rather unex- 
pected features. We will return to the discussion of this 
issue in the coming publication p2| . 
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